Diffusion at the random matrix hard edge 
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Abstract 

o 

We show that the limiting minimal eigenvalue distributions for a natural general- 
ization of Gaussian sample-covariance structures ("beta ensembles") are described by 
the spectrum of a random diffusion generator. This generator may be mapped onto the 
"Stochastic Bessel Operator," introduced and studied by A. Edelman and B. Sutton 
in [B] where the corresponding convergence was first conjectured. Here, by a Riccati 
transformation, we also obtain a second diffusion description of the limiting eigenvalues 
in terms of hitting laws. All this pertains to the so-called hard edge of random matrix 
^-j- \ theory and sits in complement to the recent work [15] of the authors and B. Virag on 

the general beta random matrix soft edge. In fact, the diffusion descriptions found on 
both sides are used below to prove there exists a transition between the soft and hard 
edge laws at all values of beta. 



m 
o . 

1 Introduction 



The origins of random matrix theory can be traced to the introduction of Wishart's en- 
sembles, matrices of the form XX' with rectangular X comprised entirely of independent 
real or complex Gaussians of mean zero and mean-square one. The spectrum of these ob- 
jects are of fundamental importance in mathematical statistics (see the comprehensive text 
P3]), and continue to generate wide interest due to their relevance to such disparate areas 
as information theory [2D], numerical analysis [7], and, along with their quaternion-entried 
counterparts, theoretical physics [27] . 

Here we consider scaling limits for Wishart-type eigenvalues at the hard edge. To explain, 
let X be n x m. If m ~ n as n j oo the minimal eigenvalues of the (non-neg ative) XX* 
will feel the "hard" constraint at the origin, while if m/n is strictly larger than one in the 
large dimensional limit, the minimal eigenvalues separate form zero and one has "soft" edge 
fluctuations (on which more below). In fact, if m = n + a with fixed a as n | oo one discovers 
an interesting family of limit laws indexed by a for the bottom of the spectrum. 
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The known results at the hard edge have thus far been based on the explicit joint density 
for the Wishart eigenvalues < Ao, Ai, ■ • • , A n _i. In particular, when X is n x (n + a) with 
integer a > — 1, that density is 
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Zf3,a 



n-l 
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fc=0 



Ko+i)-i -g: 



1.1) 



with normalizer Zp A < oo and /3 = 1, 2, or 4 for real, complex, or quaternion Gaussian 
entries. More importantly, with these choices of (3 all finite dimensional correlation functions 
of the eigenvalues are computable in terms of Laguerre polynomials (thus the common tag 
"Laguerre ensembles"). At /3 = 2 and all valid a, [22J proves the limiting distribution of the 
minimal eigenvalue is described by the Fredholm determinant of a kernel operator given in 
terms of Bessel functions, and based on this derives a second description of the limit law as a 
functional of the fifth Painleve transcendent. Other work at /3 — 1, 2, 4 hard edge include [I], 
[TO] , and • These again rely on the underlying orthogonal polynomial structure (the first 
and third reference use Riemann-Hilbert methods to replace the exponential weight e - ^/ 2 ^ 
in (11. ip with a more general e~ v ( x ' potential), and describe the eventual limit law through 
Fredholm determinants or Fredholm pfaffians. 

While the distribution on n points Ai, . . . , A n G M + defined by ( II. ip makes sense for all 
(3 > and a > — 1, the orthogonal polynomial approach breaks down outside the standard 
triple in (3. For some special choices of the parameters beyond (3 = 1,2,4, [9] was able to 
exploit the niceties of the exponential weight to obtain limit laws in terms of hypergeometric 
funcions. Still, even the existence of the general f3 hard edge limit law remained open until 
now. 

Our approach, similar to that in [TJ] , rests on the existence of tridiagonal matrix models 
for all (3. Set for any a > — 1 and (3 > 0, 



J /3,a 



VP 



X(a+n)f3 X(n-l)P 

X(a+n-l)P X(n-2)P 



X(a+2)P Xf3 

X{a+l)P 



1.2) 



in which each Xr that appears is an independent \ random variable of the indicated index. 
(We suppress here the dimension parameter n on the n x n random matrix Lp ja ). Then, 



as discovered by Dumitriu and Edelman [5], the eigenvalues of Lp >a Lp a have (11.11) as their 
joint density function. Note that, when /3 = 1 or 2, the bidiagonal (11.21) may be arrived at 
by performing Householder transformations on the corresponding "full" Wishart ensemble; 
this fact was used previously in a random matrix context by Silverstein [16J. 
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Viewing the n f oo limit as giving rise to a continuum approximation to the discrete 
operators Lp >a , an entry-wise expansion in the random \ variables led Edelman and Sutton 
to the following conjecture for the full (3 > hard edge. 

Conjecture (Edelman- Sutton [6]) Let Sk denote the k-th smallest singular value of the bidiag- 
onal operator Lp a . Then, as n j oo the family {y/nsk} converges in law to the corresponding 
singular values of 

d a 1 
dx 2yfx a//T 

in which x \— > b(x) is a Brownian motion. Here, Cp ja is understood to act on functions 
f E L 2 [0, 1] subject to /(l) = and (C p , a f)(0) = 0. 

The random Cp >a was tagged the Stochastic Bessel Operator in j6] on account of the 
zero-noise (/3 = oo) version having singular-values at the roots of J a (the Bessel function of 
the first kind). 

Our main result establishes this conjecture, though we prefer to phrase matters in a 
different way, back in terms of eigenvalues of the symmetric ensembles L^ a L^ a , henceforth 
referred to as the (j3, a)-Laguerre ensembles. Toward this, introduce the random operator of 
second order, 

<Sp ta = - exp [(a+ l)x + ^6(x)]^-|exp [-ax - ^ & ( x )]^;}> ( L3 ) 

where again b(x) is a Brownian motion and a > —1, (3 > 0. Formal manipulations will take 
you from Cp a Cf) a to <S/3, a , but the latter is better understood upon recognizing, in the spirit 
of the title, that — <25g, a generates the diffusion with (random) speed and scale measures 

2 2 
, \ -{a+l)x jsK x ) t j / j \ ax+^bix) 

■m(dx) = e vP ax and s(dx) = e vp dx. 

This motion may be built pathwise in the classical mode (see for example p2]), placing (11.31) 
on firm ground. The limiting spectral problem will require consideration of C5^ a acting on 
the positive half-line with Dirichlet conditions at the origin, and this carries over into killing 
the underlying process when reaching that point. 

Even more convenient, we may define eigenvalues/eigenvectors through the resolvent 
equation. That is, if we at first take the equation = W t° mean ip = \<8g a il>, the 

speed and scale construction provides the explicit form of the inverse, 

/oo / pxAy \ 
U 8 (dz))1>{y)m(dy). (1.4) 

Now is plainly non-negative symmetric in L 2 [lR + ,m] and the Dirichlet condition at the 
origin is automatic for solutions of ip — \Q5p\ip. Lying slightly deeper, we will see that 
(almost surely) <25o a maps L 2 [M + ,m] into C 3 / 2 " and is in fact of trace class. We have: 
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Theorem 1. With probability one, when restricted to the positive half-line with Dirichlet 
conditions at the origin, ^>p^ a has discrete spectrum comprised of simple eigenvalues < 
Aq(/3, a) < Ai(/3, a) < • ■ ■ \ oo. Moreover, with now < A < Ai < ■ • ■ < A n the ordered 
(j3,a)-Laguerre eigenvalues, 

{n\ , nXi, . . . , n\ k } {A ((3, a), Ai(/3, a), . . . , A k (/3, a)} 

(jointly in law) for any fixed k < oo as n f oo. 

Remark. The Dirichlet condition for 0^ at x = may be mapped to that at x = 1 for £^ jtl 
in the Edelman-Sutton conjecture. On the other hand, the process generated by — &p, a nas 
a natural (or free) boundary at x — +oo, which carries certain advantages over the specified 
condition for Cp >a at x — 1 in the conjecture. 

As a bit of amplification, differentiating with abandon one is led to 



along with the idea that the corresponding motion is just a Brownian motion with (shifted) 
white noise drift. In fact, modulo the multiplicative factor e x which affects a change of time, 
this is precisely the random diffusion introduced by Brox as a continuum analogue of Sinai's 
walk [3J. Theorem [1] then draws a concrete connection between random matrix theory and 
the lifetime of this random process in a random environment which has been the subject of 
continued investigation since its introduction (see fTS] and the many references within, or 
the recent [2] for a spectral point of view). 

Our second description of the hard-edge is a corollary of the first, employing Riccati's 
map to transform a solution of, the suitably interpreted, <5/3 ja i/j(x, A) = \i[)(x, A) for any fixed 
A > into one of 



understood in the sense of Ito. The point is: Sturm's oscillation theorem implies that the 
eigenvalues of (5^ are counted by the zeros of A), and those zeros correspond to places 
where p = ip' /ip, which solves (11.51) . hits — oo. 

Theorem 2. Let P x c denote the law of p(-) = p(-\ a, j3, A) starting from position c at time 
x. Let also v x {dc) = P Xj+00 (m G dc) where m is the passage time of p to — oo. Then, 




dp(x) = -jjjp(x)db(x) + ( (a + V)p{x) — p 2 (x) — Ae x ) dx, 



(1.5) 



P{A ((3,a)>\) 
P(A k (P,a)<\) 




f ({oo}) and, more generally, 
/ v (dx 1 )v Xl {dx 2 ) ■ ■ -v Xk {dx k+1 ). 
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Theorem 1 is proved in Section 2, and Theorem 2 in Section 3. We conclude the intro- 
duction by describing a general transition between the hard edge laws just described and the 
form of the (3 > soft edge laws established in [15] . 

Remark. It is natural to ask whether from Theorems 1 and 2 one might recover the Painleve 
or Hypergeometric descriptions of the hard edge from [22] or [9] respectively, and then go 
further by finding explicit formulas of the distributions at all (3 > 0. Thus far the answer 
is no, even when (a + 1) = 2/(3 and A (a, (3) is just an exponential random variable (as 
is easily seen from the joint density fll.ip ). With that choice of parameters, the speed 
m'(x) = e ^ a+1 ^ x VP hx ^ f (g ^-diffusion turns out to be a martingale, but we do not see 
how to make use of this. 

Soft edge and transition 

The random matrix soft edge corresponds to the scaling limits of the maximal, rather than 
minimal, eigenvalues in the Laguerre ensembles. Historically, these laws were discovered 
first by Tracy and Widom ([21] and [23]) in the context of a different class of random 
matrices, the Gaussian Orthogonal, Unitary, and Symplectic ensembles. The latter are nxn 
real symmetric, complex hermitian, or quaternion self-dual matrices with Gaussian entries. 
There is again an explicit joint spectral density, which takes the form of a constant multiple 
of Yl ~ A J '| /3 e~( /3 / 4 )^ A ? with /3 = 1,2 or 4 respectively. And once more, all correlations 
are given in terms of orthogonal polynomials (now Hermites). Tracy and Widom proved 
that the appropriately scaled largest eigenvalues have distribution functions described by 
Painleve II (via a more basic formulation in terms of Fredholm determinants/pfaffians of an 
Airy kernel). Matters were later carried over to the Laguerre soft edge by a collection of 
authors. 

As before, one may consider the general "/3-Hermite" laws. [5] provides a separate family 
of tridiagonal matrix models for these laws (though see also [21] for an earlier application 
at (3 = 1,2 ), and in direct analogy with Theorems 1 and 2 the authors and B. Virag have 
previously proved: 

Theorem [Theorems 1.1 and 1.2 of [15]) The largest eigenvalues in either the (3-Laguerre 
or (3-Hermite ensembles have scaling limits given by the law of the top eigenvalues of the 
random Schrddinger operator —Tip = J-^ — x + -^=b'(x). There is also a description of the 
limiting soft- edge eigenvalues equivalent to that in Theorem 2, with the diffusion 

dp(x) = -^db{x) + (A + x -p 2 (x))dx (1.6) 

in place of (11.51) . 
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The first part of this result, the identification of —Tip at the soft edge, proves a different 
conjecture of Edleman and Sutton from [6j. For obvious reasons, we refer to the distribution 
of the top eigenvalue of — Hp as the general beta Tracy- Widom law, notated TWp. 

Returning to the (/3, a)-Laguerre ensembles, it is well understood that if a tends to infinity 
with n so that \mx n -* 0O m/n = lim n ^ 00 (ra + a)/n > 1, the limiting spectral measure is pulled 
away from the origin and one sees soft-edge behavior at both the minimal and maximal 
eigenvalues. Thus, one expects that by taking a — * oo after n —>■ oo, the hard-edge becomes 
a soft-edge and creates a link between these families of distributions arising in random matrix 
theory. Borodin and Forrester pQ have shown that this is indeed the case in the classical 
(3 — 1, 2 and 4 settings. Their work rests on the aforementioned determinantal forms of the 
underlying distribution functions. Employing just the diffusions ( 11. 5ft and (11. 6ft . our final 
result shows that this transition holds true at all /3 > 0. 

Theorem 3. With A (/3, a) the limiting smallest eigenvalue in the a) -ensemble and TWp 
the general beta Tracy- Widom law, 

ft l 

as rj — > oo. 

Theorem [3] is proved in Section 4. 

In summary, together with [15j the present provides a complete picture of the extremal 
laws of random matrix theory, at all values of the natural parameters. This leaves apart 
the general f3 spectral bulk, which has recently been treated by Valko-Virag [25J (for the 
/9-Hermite ensembles) and Killip-Stoiciu [13] (for the circular (3 ensembles, generalizing the 
eigenvalue laws for the Haar distributed unitary group). 

2 Convergence of the spectrum 

The key is to prove the almost sure strong convergence of the resolvent operators, or really 
a similarity transformation of the sequence of (L ( g a Lj (l )~ 1 matrices to a version of &g a - As 
we will see, all these objects may be viewed as integral operators with well-behaved kernels, 
allowing for an efficient verification of the necessary compactness. 

Outline 

Set M/3 >a = SLp^S^ 1 where S is the anti-diagonal matrix of alternating signs Sij = (— l) l 5j + j_ n _i. 
The spectrum is unchanged (we may work with Mp^MT rather than L^LT), and we 
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record 



M, 



f3,a 



1 

7? 



X(a+l)f3 
~XP X(a+2)/3 

~X2/3 X(a+3)/3 



—X(n-l)P X(a+n)l3 

where the additional notation is intended to emphasize the independence of the processes 
along the main and lower diagonals. 

Wishing to track inverses, we first note the readily checked fact: 

Lemma 4. For any lower bidiagonal matrix B = b it j (that is, bij = if j > i or j < i — 1), 

the inverse, when it exists, is lower triangular and has the expression 



1,3 



[-l)^ 1 T-r b k+ i k . 

[[— for j < i. 



hi bk,k 

k=j 



Next, observe that for any A = a i: j G 
L 2 [0, 1] which does change the spectrum: 



(Af)(x) 



= 2^ Oijn \ 

7=1 Jx 3-1 



there is a natural operator embedding into 



f(x)dx for < x < Xi, 



where hereafter we define Xi — i/n, for i — 1,2, ... ,n. Thus, moving attention to (nMp^Ml ^ ' 



Pa' 



(after introducing the appropriate hard-edge scaling), the action of n 1//2 Mg* on L 2 [0, 1] 
reads 



\nx\ j-t, — \nx J — 1 _ 

((v^*0-7)(*) = £d^: II u 



f(x)dx. 



(2.1) 



J= l X([nx\+a)f3 X{k+a)f3 J x^ x 

In other words, n~ x l 2 Mp \ is equated with the integral operator K^ a with (discrete) kernel 

'pn 



k p,a( x >y) = exp <^ VlogXfc/g - logX(fc+a)/3 > 1l(x, y) 

X(i + a)f3 J 

where 1 L = l{ Xl _ 1 <x<x l }^{x^ 1 < y <x J } and % > j. 

With this set-up, the basic convergence result we need is the following. 

Lemma 5. There is a Brownian motion b(-) such that for x < y lying in (0, 1] 



X(\nx\+a)l3 Vx 



and 



[nx\ 



( lo §^/3 ~ lo gX(fc+a)/3) => (a/2) \og(y/x) + 

k= \ny\ 



db z 
tl VP?'' 



(2.2) 



(2.3) 



in law in the Skorohod topology. Morever, the exist tight random constants n n > and 
k' > which are independent of j3 so that 



sup < n n (2.4) 

l<fc<n X{k+a)(3 



and, with T{x) = 4 log 



i-1 



(\ogXkp ~ logx (fc+a)/3 ) - (a/2) log(jA) < <(1 + T 3 /\ Xl ) + T 3 /\ Xj )) (2.5) 

k=j 

for all 1 < j < i < n. 

The first part of the lemma ( fl2.2[) and (12.31) together) identifies the limiting operator 
Kp a . Namely, for n f oo it should be that kp a (x, y) approaches 



kf3 ja (x,y) = x 2%xp 



db* 



y a/2 l y<x . (2.6) 



The second part, or the bounds (12 ,4p and (12.51) . provide the needed compactness and more. 
As we shall prove: 

Lemma 6. Kp A is almost surely Hilbert- Schmidt. Also, there exists a probability space 
on which all K^ a and Kp A are defined, and such that any sequence of the operators K^ a 
contains a subsequence which converges to Kp^ a in Hilbert- Schmidt norm with probability one. 
In particular, for whatever n f oo we can find an n' f oo along which 

lim / / \kp {x,y)(uj) - k^ a {x,y){u)\ 2 dx dy = 

" Too J Q J Q 

almost surely. 

Granted this we may complete the proof of the main result. 

Proof of Theorem [IJ Working on the probability space promised in Lemma [61 the argument 
is reduced to a deterministic setting. Start with the scaled minimal a)-Laguerre eigenvalue 

nAo(n) = inf (v,nM^ a Ml a v) = ( sup (f^K^fK^J))' 1 = \\(K% >a ) T K^ 1 , 

where || ■ || is the L? \— > L? operator norm, and the final equality holds simply because 
(Kg a ) Kg a is non-negative symmetric. Assume for the moment that, as claimed, Kp a is 
almost surely Hilbert-Schmidt. Then Kp a K/3 ia is non-negative symmetric and compact (trace 
class even) with a well defined maximal eigenvalue also equal to the norm ||i^J a i^ j0 ||. For 
short, we notate A > Ai > ■ • • the eigenvalues of Kj a K/3 >a = B, and similarly write {A£} for 
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the (decreasing) eigenvalues of (Kg a ) T Kg a = B n . The simplicity of the limiting eigenvalues 
is also assumed here; it will be deduced from the differential form of the eigenvalue problem 
introduced in the next Section, see 13.11 and the surrounding discussion. 

Next, for whatever sequence n f oo, Lemma [6] allows a choice of subsequence along which 
| \Kp' a — Kp >a \ \hs 0. The same holds for the transposes, and hence B n i converges strongly 
to B. It follows that the norms themselves converge along this subsequence: Aq = ||P n '|| — > 
\\B\\ = Aq with probability one. But this is to say that for any sequence of the original 
eigenvalues n\o(n), there exists a subsequence along which these points converge almost 
surely to 1/A . That is of course equivalent to the full convergence statement. 

As to nAi,nA2, we first show that the convergence (along perhaps a further subse- 
quence) of the ground state eigenvectors is a by-product of the above. Define {f n } and / in 
L 2 [0, 1] with unit norm by 

(f n ,B n f n )=A n Q , (f,Bf)=A . 

Remaining in the introduced setting, we have ||-B n /n||L 2 — ► H-B/Hl 2 - Also, being uniformly 
bounded in L 2 , f n has a weakly convergent subsequence: f n i — ^ /oo. Then, for any G L 2 , 

(0, B n f n - BU = (0, {B n - B)f n ) + (50, f n - U 

tends to zero (the first term by norm convergence, the second by boundedness of B) . Having 
weak convergence (B n f n — ^ Bf^) plus convergence of its norm, we conclude there is a 
strongly convergent subsequence of {B n f n }. Coupled with B n f n = Aq/„ and Aq — > A , this 
implies a strongly convergent subsequence for the {f n } themselves, which by continuity can 
only wind up at /. 

Finally, place yourself along this sequence where | \f n — f\ — > 0, and denote by P/ n the 
projection onto the orthogonal complement of f n in L 2 . At once we find that Pf n B n Pf n con- 
verges strongly to PfBPf (with obvious notation), and so A" = \ \Pf n B n Pf n \ \ — > ||Pj5P/|| = 
Ai. The implication for n\\ is clear, and an induction argument extends the picture to the 
almost sure convergence of any finite number of Laguerre eigenvalues. 

Reflect upon the fact that we have proved convergence (in law) of say {(^A^) -1 } for 
k = 0, . . . , m to the top m eigenvalues of the integral operator B = K^ a Kp a which we now 
write out. In particular, its spectral problem reads 

f(x) = X [ x a/2 e S *^y- {a+1) [ l ' e ! *^z a/2 f{z)dzdy (2.7) 

Jx J0 

= \f (xy) a/2 ( I e~ 2S *vm z- {a+1) dz\e f *^e f y^f(y)dy, 

JO \JxVy J 

after an integration by parts. Again, we seek here an / 6 L 2 [0, 1], which inherits the 
continuity and vanishing of the kernel at x — 1. 
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ri dbg 

To recover the advertised limit operator (II. 4p . make the substitution g(x) = x~ a ' 2 e Jx Jp* /( 
in conjunction with the time change f 1 s~ l l 2 db s = fo(log(l/a;)) with a new Brownian motion 
b to express (12.71) in the equivalent way: 



i / rl 



g {x) = A / / e-^ b{logl/z) z-^ dz ) g(y) y a e^ h ^ 1/y) dy. 

J0 \JxWy 



With / e L 2 [0, 1], g resides in L 2 ([0, l],m) for m{dx) = x a e^ b ^ logx ) dx. Last, the change 
of variables (x,y) i— > (e~ x ,e~ y ) will produce the form of (11. 4p quite exactly, along with 
transforming the Dirichlet condition at one (/(l) = <?(1) = 0) into that at the origin (?/>(0) = 
(poexp)(0)=0). □ 

Estimates 

Before establishing Lemmas [5] and [6] we make the simple observation: 

Proposition 7. For any constant C and a > —1, the integral operator on L 2 [0, 1] with kernel 
k c {x,y) = Cex V [c(log(l/x)) 3 / 4 + C(\og(l/y)f 4 



.(a+l)/2 1 ^ <x 



is Hilbert- Schmidt. 

Proof. The change of variables x = e~ s and y = e~ l employed just above produces 

•1 pi POO poo 



k c {x, y)\ 2 dxdy = C 2 / e 2Cs3/4+as / e 20 ^-^ 1 dt ds 



'0 JO JO 

and the latter is clearly finite if (and only if) a > — 1. □ 

Proof of Lemma\^ Now taking Lemma [5] for granted, we can find a subsequence over which 
we have the joint convergence in law, 

VnP ( 1 



X([nx\+a)l3 \V X 

\nx 



, < x < 1 , 



/ f x db \ 

(logXfc/3 - logX(fe+a)/3) ( (a/2) log(i//a;) + / < y < x < lj , (2.1 

= i IT?/ 1 y * 



k=\ny\ 

I / 

At„ , ft„ =7 > " Hi. Hi . 



Then, Skorohod's representation theorem (Theorem 1.8, Chapter 2 of [8]) furnishes a prob- 
ability space on which each of the above occurs with probability one. The first two items of 
(12. 8j) take place a.e. in (0, 1], and so on this new space it holds 

P(\imkl a (x,y)(uj) = kf,, a (x,y)(u) for a.e.x,y E [0, l] 2 ) = 1. (2.9) 

Vnjoo / 
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That 

-1 /■! 



JO 



\tf3,a( x iV)( u ) - kj3, a (x,y)(u)\ 2 dxdy -> a.s. 
will follow if we can supply an a.s. finite constant C(oj) such 

sup k% a (x,y)(uj) < kc{u)(x,y) and kp, a (x, y)(u) < k c(uj) (x, y) (2.10) 



n>0 



for almost all x, y G [0, 1] and cj. 

Again by Lemma El for each n it holds 

v)(*«) < M^K~ (a+1) V 2 ex P K(^)(! + t3/4 (^) + t3/4 (%0)] (2.ii) 

where x G [xj.x i+ i),?/ G [t/j, 2A/+i) - But now we are allowed to assume that both K n and n' n 
converge, and thus are bounded almost surely by say 2(n V k') for sufficiently large n. The 
continuity of the functions x \— > T(x) and x i— > x p (on (0, 1]) then enables us to fit the right 
hand side of ( 12.1 ip under a fixed kc independently of n. For the limit kernel, kp ja (x,y)(uj) 
simply note that its exponent could have been expressed from the start as 

" x db, 1 



6(log(l/x))-6(log(l/y)) 



(The equality is in law with a different Brownian motion living on the same probability 
space.) By the law of the iterated logarithm, b(a) < c(u)(l + [a log log (1 + a)] 1 ^ 2 ) for a 
random c(u) and all a > 0, and certainly [a log log (1 + a)] 1 ^ 2 < c'a 3//4 with a (non-random) 
d and all a large enough. Thus, the second half of (I2.10p holds with C(u>) < j3~ l l 2 c' c(u), 
and the proof is complete. 

Note here we immediately passed to a fixed subsequence and then chose a favorable 
probability space, while the statement of the lemma was worded with the convergence of the 
(P, a)-Laguerre eigenvalues in mind. That is, build all kp a , each tied to a (/?, a)-Laguerre 
eigenvalue, on the same space as kp >a and then note for whatever n f oo there is a subsequence 
along which everything above holds. Either way the upshot is the same. □ 

Turning to the proof of Lemma El we record without proof the following facts. 

Proposition 8. For Xr a chi random variable of index r > 0, 

E[xl] = 2P nM (2 - 12) 

for any p > —r. Also, as r — » oo ; 

E[\og X r] = ilogr-|- + 0(l/r 2 ), Var[\og X r] = ^ + 0(l/r 2 ), (2.13) 
while E[(\ogXr — -ElogXr) 2m ] = 0{l/r m ) for positive integer m. 
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Proposition 9 (After Theorem 1.3, Chapter 7 of |8J). Let y n ^ be a sequence of mean- zero 
processes starting at with independent increments Ay n ^- Assume, 

nE(Ay n<k ) 2 = f{k/n) + o(l), nE{Ay n . k f = o(l) (2.14) 

uniformly for k/n in compact sets of [0, T) with a continuous f G Lj oc [0,T). Then y n (t) = 
Un,[nt] Jq f 1 ^ 2 (s)db(s) with a standard Brownian motion b (in the Skorohod topology). 

Proof of Lemma{3 Start with (12.31) . By the first estimate of (12.131) . 

\nx\ 

tim ( Elo SXk/3 ~ E\ogX(k+a)f3) = 2 l °g(y/ x ) 

k= [ny\ 

uniformly for y < x restricted to compact sets of (0, 1]. Thus, for ( 12. 31) it is enough to 
demonstrate the weak convergence 

n „1 

(log X(k+c)p- Slog X(k+c)p)^ / (2/^)" 1/2 db(z) (2.15) 

k=[nx] Jx 

where c is any fixed number. Indeed, the exponent of the discrete kernel is comprised of 
two such independent sums, and the promised limit will follow as &i + 62 = "\/2&3 in law for 
independent Brownian motions b%, 62, &3- Now refer to Proposition [9] and view the processes 
on the left of (12.151) as starting from at x — 1 and evolving toward x = (or take t — 1 — x 
in the proposition). Then, the second estimate of (12.131) yields the first part of (I2.14p with 
f(t) = l/(2/?t); the estimate right after (12.131) with m = 2 produces the second half of (12.141) 
clS X IS always > 0. This finishes the job. 

The convergence (12.21) is easier. For any fixed x G (0, 1], it is just an instance of the law of 
large numbers. The tightness required to ensure process level convergence is also elementary: 
via (I2.12p one can obtain the increment bound 





/ r(3 \ 


1 X( r+a+ l)f3 V 


X( r +a)f3 ' 



which more than suffices. While here we dispense of (12.41) . First use the sum bound, 

p ( sup > m) < jrp ( < 1 



l<fc<n X(k+a)/3 J ~^ \ VW M 

Then, employing the explicit density P(xr £ ds) = ^72) s r ~ l e~ s2 / 2 ds, one can perform a 
Laplace-type estimate to find the k-th term on the right hand side is upper bounded by 
C(^/~e/M) k with C depending only on (3. Since Y^k=i(Vz/M) k ma y t> e m ade arbitrarily 
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small by choice of M, the desired tightness of the random variables sup k< n (V^P I 'X(k+a)p) 
follows. 

The final piece, or (12.51) . is the most elaborate but really comes down to reworking the 
standard proof of the upper bound in the law of the iterated logarithm. Define, 

71-1 

A x = ^20°SXk/3 ~ !ogX(fc+a)/3) - 2 logO'/n) 

k=j 

for x e [xj,Xj + i), and h(x) = [2x log logx] 1 / 2 . We will in fact show that 

sup ((A™. V 0)/h(T(xj))) are tight in distribution, (2.16) 

where again T(x) = 4 log ~. This is stronger than what is claimed. 
Set 



n\Ai-l 



v /£- + l\ a/ 

i- s «p(^) = n^- ^ . ^d z» s {y^e\{y; 

k=j X(k+a)/3 \ K J 

with a small positive A (the precise conditions on A follow shortly). The sequence j i— ► Z"_j 
is a martingale for j = 1,2,... with E[Z^\ = 1 for all j. Hence, by Doob's inequality 

P( max > e Xb ) < e~ x \ 

\£<j<n-l J J 



or 



P( max (XAl - logS[exp(A^)]) > b) < e^ b (2.17) 

\t<1<n— 1 J J / 



for 6 > 0. 

For the next move we need an estimate on the moment generating functions of A™., the 
proof of which we will return to at the end of the section. 

Claim 10. For all A > sufficiently small (A < (/?/2)[(l + a) A 1] will do), 



E 



e i 



expj^Ml/x.O + e^j)} (2.18) 

with |6 n (i)| < C for constant C = C(a, j3). 
Using fl2TT8D in fl2TT7jl . we have 

p(wp^ - Yp Xog{l,t) + x 0n(nt) } - h,x ) - e ~ Xb 

with n (i) understood via interpolation. Now choose > 1, a positive constant M and set 
A = M9~ m h(9 m ), b = Mh(0 m )/2. (To choose M large one must take B large as well to 
respect the condition on A set down in Claim [TOl ) The previous display will then imply 

sup A n t > (M+ l)h(6 m ) ) < (m\og6)~ M2 . (2.19) 

i <T(t)<6» m + 1 
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Here we have used the uniform bound on Q n (t) to fit A 1 Q n (t) under h(8 m ) by choice of 9 
and so M. In particular, A" 1 = M- l 6 m h- 1 (6 m ) < M~ l h(6 m ) if log log > 1/2. 
Finally return to the goal (j2.16j) . re-expressed as seeking a bound of type 

P ( sup [A?] + /h{T(t)) > N j < e(N) where e(N) j as N | oo. 
\o<t<i / 

Note in addition that the supremum inside the probability over any truncated range x < t < 1 

(for x > 0, rather than < t < 1) poses no problem. Indeed, the process 1 1— > A" has already 

been shown to be convergent in that regime. On the other hand, the troublesome tail is 

bounded by 



m=l 



sup A?/h(T(t)) > N 
i <T(t)<e m + 1 



< y^(mlog^) 

m=l 



-{N-l)' 2 



with the aid of (I2.19p . completing the proof. 

Proof of ClaimM By (12121) . the left hand side of (12151] equals 



□ 



n-l 

II r 

fc=j 1 V 2 y 



p / (fc+q)/3 



Hl\ Aa/2 



fe 



Taking logarithms, we must estimate the sum ^2^=) s k where 



k=j 



s k 



io g r 



kp + A 



logT 



(fe + o)/3 - A 



-iogr 



(fe + a)f3 



Xa 
+ ylog 



(2.20) 



Introduce Stirling's approximation in the form 



logT (z) 



log z + z — 



log 2tt 



< 



c 

^2' 



2 12z 

The 0(,z~ 2 ) error term produces a constant multiple of fe -2 when applied in (12.201) . Dif- 
ferences such as (fe/3/2)" 1 — ((k/3 — A)/2)~ : and the like stemming from the 1/(122) terms 
are similarly bounded. When summed, both contributions produce constants which are 
then absorbed into the n . Also, the constant and z-terms obviously cancel throughout the 
log-gamma expressions when the above estimate is applied in (I2.20p . 

Move to the terms of type (z — 1/2) logz. A bit of algebra will lead to 



k/3 



log 



a(3 
+ T log 



1 + 

A 



(fe + a)/3 



A 



(fe + a)(3 
2 l0g ( 1+ fe 



A. 

2 log 

Xa 



1 + 



-) 

k/3 J 



/ I 



A 



(fe + a)l3 



— log[ 1 + - ] +0[- 



Since | log(l + s) — s\ < s 2 for s > —1/2, we conclude = ^ + 0{k 2 ), which establishes 
the claim upon summation from j to n — 1. □ 
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3 Riccati map and a second diffusion 



Riccati's substitution takes a linear second order operator into one of first order, at the price 
of introducing a quadratic nonlinearity. Its use in the study of random spectra has a long 
history, dating back to Halperin [11]. To employ it here we must first recover the differential 
form of the eigenvalue problem from the established integrated version ip = A(5^^, which 
reads in full: 

poo pxAy 

ip(x) = X J J s{dz) ip(y) m(dy) 
= X 



JO 

oo / pxAy 



exp [az + 4%b(z)) dz ip(y) exp [-(a + l)y - 4*b(y)] dy. 
'o \Jo / 

Noting that any / e L 2 [m] is also in L^m], &p a f is easily seen to be differentiate after 
writing the right hand side as separate terms. This property is inherited by ip, and we 
compute 

poo 

ij/(x) = Aexpfax + ^b(x)] / if>(y) exp[-(a + l)y - ^b(y)] dy, 

J X 

to find that ift is actually in C 3 ^ 2 ~ . Continue by taking (Ito) differentials to arrive at the 
system 



dif)'(x) = -^ip'(x)db(x) + (J a + %)if)'(x) - Xe x ip(x)j dx, 
dip(x) = ip'(x)dx, (3.1) 

which is the appropriate way to interpret tSp^ip — Taken independently of the preceding 
developments, (13.11) has globally Lipschitz coefficients of linear growth, and as such defines 
(for fixed A) a unique Markov process x i— »■ ijp{x), ip'(x)) for any specified (^(O), ^'(O)) pair. 
As the uniqueness is pathwise we conclude along the way the simplicity of the corresponding 
eigenvalues: for given A, any two L 2 solutions of ip = AG5^?/> vanishing at the origin must 
be constant multiples of one another. 

Now bring in Riccati's map, p(x) = ip'(x)/ip(x), valid away from the zeros of ip. Since ip 
is continuously differentiable, we find from (13. 1 p and elementary calculus: 

dp(x) = -^p(x) db(x) + ({a + f)p{x) - p 2 (x) - Xe~ x ^j dx, (3.2) 

defining yet another Markov process for any fixed A. The relevance of (13.21) in counting 
eigenvalues of (S^a is first understood through the truncated operator £5^, indicating (25^ 
restricted to [0, L] with Dirichlet conditions at both endpoints. 



Lemma 11. Consider the unique diffusion p(x) = p(x; A) started at +oo at x = 0, and 

*> L 



restarted at +oo immediately after any passage to — oo. The number of eigenvalues of <25f 



less than X is equal in law to the number of explosions of p before x = L. 
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Proof. This is well understood, and our treatment here is much the same as in |15j . Section 3. 
Take the sine-like solution of (13. II) . that is, i/)q(x, A) subject to ipo(0, A) = and ip' (0, A) = 1. 
Plainly, A is an eigenvalue of <8p a if only if ip (L,A) = 0. Regarding the ground state 
eigenvalue A (L): if for any A i/)o(x, A) > for < x < L, then it must be that A (L) > A, 
as an examination of (13.11) shows. That is, the event that {Aq(L) > A} is equal in law to 
the event {x i— > ipo(x,X) has no roots before x = L}. Continuing, additional zeros of the 
(almost surely continuous) function A i— > i/)q(L, A) (and so additional eigenvalues) only occur 
by increasing A, whereupon all other roots (in the x-variable) move to the left. This equates 
the event that the k-th eigenvalue of (J5^ a lies above a fixed A and the event that i[>o(x, A) 
has at most k — 1 roots on (0, L). 

Now move to the p(x,\) formed from tpo(x, A) and its derivative. By appealing again 
uniqueness of solutions to (I3.ll) . note ip and ip' Q cannot vanish simultaneously. (In particular, 
the zeros of ?/>o are isolated, and must be either finite in number or form a sequence tending 
to infinity.) Thus, at any root m of x i— > ipo(x, A), including m = 0, an examination of signs 
shows that lim e jop(tn + e, A) = +oo and, when m > 0, lim £ j p(m- — £, A) = — oo. That 
is, counting roots of i/jq(-, A) is to count passages of the corresponding p(-, A) to — oo, after 
subsequent re-starts at +oo. 

To see that the p-picture stands on its own is to show that there is a unique solution of 
(13.21) starting from +oo. Replacing the — Xe~ x term in the drift with any negative constant 
produces a homogeneous motion with an entrance boundary at +oo (and which hits — oo 
with probability one). This process (begun at +oo) may be constructed unambiguously via 
speed and scale, see again [12J. By successive dominations of the inhomogeneous p in the 
statement by such homogeneous versions over all short times, one may conclude the existence 
and uniqueness of the former. □ 

Theorem [2] now follows by taking L — » oo in Lemma [TT1 with the aid of the next fact. 

Lemma 12. As L — » oo ; the top k eigenvalues of<3p a converge to the top k eigenvalues of 
(3/3^ a with probability one. 

Proof. This again demonstrates the advantage of having explicit inverses. Now (£5^ J" 1 acts 
on L 2 ([0, L], m) via 



Plainly, sl(x, y) < JJf v s(dz) and lim^oo sl(x, y) = J* v s(dz) pointwise in x and y, almost 




where 




xAy 



fxvy < dZ ) 

Jo s ( dz ) 



{x,ye[0,L]}- 
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surely. By dominated convergence we have in the same mode that 

oo poo poo poo / pxAy \ 

/ f{x)s L (x,y)g(y)m{dx)m(dy) -> / / f(x) ( / s(dz)) g(y)m(dx)m(dy) 
o Jo Jo Jo \Jo J 

for all /, g G L 2 [IR + ,m], and 

tr((S^ a ) _1 = / s L (x,x)m(dx) -> / / s(dy)m(dx) = trtS"^. 
Jo Jo Jo 

But these last two items imply convergence of to (5^ ja in trace norm (see [17], Theorem 
2.20); the convergence of the eigenvalues then stems from the same style of argument used 
in the proof of Theorem [IJ □ 



4 The Hard-to-Soft transition 



Borodin- Forrester pQ discovered a transition between the hard and soft edge distributions at 
(3 = 1,2, and 4. Their proof rests on the explicit Fredholm determinant or Fredholm pfaffian 
form of these laws. For example, at (5 = 2 one has that 



P(A (2,a)> \) = 1 + J2 



k=l 





f x 


/ dx\ ■ ■ 


■ / dxk det 


Jo 


Jo 



K Besseli-Eii % 



- i,j=l,...,k 



(4.1) 



while 



P{TW 2 < A) = 1 + 



fe=i 



POO 


POO 






/ dx\ ■ ■ 


■ i dxk det 















(4.2) 



Here, 



K 



Bessel 



x-y 



with J a the usual Bessel function of the first kind, which is replaced by the Airy function in 

Ai(x)Ai'(y) - Ai'(x)M(y) 



KAiry{X-, y) - z 



x-y 



For (3 = 1 or 4 the determinants in (14. ip and ( 14.21) are replaced by quaternion determinants 
(or, equivalently, pfaffians), but are comprised of the same class of functions. Further, it is 
a fact that, suitably scaled, J a goes over into the Airy function as a —>■ oo, and the analysis 
of pQ demonstrates that one may pass this limit inside the various multiple integrals in (14.11) 
and its analogues. 

By a much different method, employing the Riccati correspondence, we show the same 
type of phenomena holds at all (3 > 0. 
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From Theorem [21 the event that {A (/3, a) > A} is equivalent in law to the process 



dp(x) = -^p(x)db{x) + ((a + ^)p(x) -p 2 (x) - Ae x j 



dx 



never hitting — oo. While from [15] we know that the probability of the event {TWp < /i} 
equals the chance that a separate motion q given by 

dq(x) = -^db(x) + (x + jj, — q 2 (x))dx (4.3) 

also never hits — oo. (Both processes are begun at +oo.) The question is then: with the 
scalings 

2 

a = — — > — 1 and A = rj — i] 2 ^ 3 fi, 

P 

does the chance of p-explosion go over into that of a g-explosion for large rp. 
To understand the mechanism, set /x = for a moment. This scaled p solves 

dp(x) = -^p(x)db(x) + (2^/rip(x) — p 2 {x) — 7]e~ x )dx, 

and obviously p — pf y/rj explodes or not with p while satisfying 

dp(x) = ^p(x)db(x) + y/rj(2p(x) - p 2 (x) - e~ x )dx. 

For rj | oo, p comes quickly to the place p — 1, and, if it manages to tunnel through this 
point in a short time, explosion is hard to avoid. Within this excursion from 1 + to 1~ in 
a small x-window, the g-motion emerges. To make this explicit we will use the following 
convergence criteria. 

Proposition 13 (After Theorem 11.1.4 of [H]). Let a(t,z) and b(t,z) be continuous from 
[0, oo) x M into R. For each w G M, let the solution of the martingale problem for a and b 
(diffusion and drift coefficients respectively) begun from w at t = s be unique. Denote this 
solution by P SjW . Suppose next that there are {a n } and {b n } satisfying 

sup sup sup (\a n (t,z) \ + \b n (t,z)\) < oo 

n>l t<T \z\<M 

and 

r 

n—>oo 



lim / sup (\a n (t, z) — a(t, z)\ + \b n (t, z) — b(t, z)\) dt = 

'0 \z\<M 



for all T > and M > 0. Then, if P™ w is a solution of the martingale problem for a n and 
b n starting from (s,w), P™ w — > P s 



s.w ■ 
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Proof of Theorem Restoring a generic value of \i we write 

dp{x) = ^p{x)db{x) + r] 1 / 2 (2p{x) - p 2 (x) - (1 - 7T 1/3 /i)e~*)cfe. (4.4) 

Here p(0) = +oo, while to utilize the proposition it is convenient to move the starting point 
to a finite place. 
Certainly, 

P+oo (p never explodes j > P\ +£ i^p never explodes j 
for whatever e > 0. Also, 

P+oo (p never explodes j < P +00 ^p never explodes, rrii +e < 5 ) + P +00 (mi +e > 5 ) 

where m c is the fist passage to the point c and 5 > 0. By the Markov property and monotonic- 
ity, the first term on the right is less than the (P^i+^-probability of no explosion. We wish to 
bound the second term from above for large rj, and to that end note that P +00 (m a < m 5 a ) = 1 
where is the passage time of the homogeneous process ps in which the appearance of e~ x 
in the p drift is replaced by e~ s . (The obvious coupling is used.) Hence, 

P+oo (mi +e > d) < -P +00 [m^ +£ ] = -j j s(dy)m(dx) (4.5) 



for m(dx) and s(dx) the speed and scale measures of p^: 

m (dx) = -^-e-^ x) dx, s(dx) = e^ {x) dx, ip(x) = - 
p x 2 2 

and c v j = (1 — iM]~ 1 ^)e~ s . Next choose 



x — 2 lnx — c„ g~ 1 

x 



£ = £( v ) = Mrp 1/6 , 5 = 8(rj) = ^?T 1/3 , (4-6) 

K 



where K > 1 and M > \J\[i\ + 2. These last precautions imply that ip(x) is increasing for 
x > 1 + e. Then an exercise in stationary phase allows the continuation of (14.51) as 



/* OO "1 px 

P + Jm 1+£{v) > Siv)) < Krj 1 / 3 / - I e -VnM*)-m] dydx < C - 

V ' Jl+Mr)- 1 ' 6 X Jl+Mv- 1 / 6 



K_ 

I+MJ7-V6 X* Ji+Mr]- 1 / 6 M 

for 7] I 00 and a constant C depending only on (3, the inner integral concentrating at the 
upper limit y = x. In summary, for p paths we have that 

Po,i+e( v ) (m-oo = 00) < P ,+oo (rn-oo = 00) < Ps( v ),i +£ ( v ) (m-oc = 00) + C— (4.7) 

holds for all large 77. 
Now bring in 

q v (x)= V 1 / e (p( V ~ 1 / 3 x)-l), 
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and note that, when p begins at (0, e(i])), q v begins at (0, M), and when p begins at 
(S(r)), e (77)), q v begins at (K~ l , M). Further, q v hits —00 if and only if p does, and a substi- 
tution in f)4.4p shows that q n satisfies the ltd equation 



dq v (x) = ^ [l + 7]- 1/6 q v (x)} db(x) + \-qfa) + V /3 { 1 - (1 - r]^ 3 ^e^ 3 



dx 



with a new Brownian Motion b(x) = r/ 1//6 6(^ 1//3 a;). Given unique strong solutions in both 
instances, Proposition [13] easily applies with 

a v (t, z) = (2//3)[l + ^z} 2 and b v (t, z) = [-z 2 + ^{l - (1 - Tj-^e^- 1 ^)], 

the ^-coefficients, and a(t, z) — 2/(3 and b(t, z) = —z 2 +fi + t, the g-coefficients (recall (14. 3p ). 
That is to say, lim^oo E x iC [4>{q v )\ = E x>c [(f>(q)] for all bounded continuous functions of the 
path, and, by approximation we also find, via (14. 6 p and (14.71) . that 

Po,a/ (q never explodes j < lim inf P ,oo (p never explodes j 

< lim sup Po.oo (p never explodes") < Pjt^^iq never explodes") + C— . 

Note while q 1— > m_ 00 (g) is not continuous, q t— > m_i(g) is for any L finite (outside a set of 
measure zero). It follows that we have the distributional convergence of m_z / (q , ?? ) to m_L(g). 
The approximation required above is then to show that: lim^oc m_L(g) = m_oo(g) holds 
in probability, with the same limit taking place uniformly in r] when q v replaces q. That 
all processes involved have exit barriers at —00 makes this routine. To finish the proof, 
let M and then K tend to infinity. The g-law is continuous in its initial time, and that 
liniM^oo P c ,m = P c ,oo is a byproduct of +00 being an entrance point. □ 
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